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ABSTRACT 



Magnetically-driven non-stationary acceleration of jets in AGNs results in the leading parts 
of the flow been accelerated to much higher Lorentz factors than in the case of steady state accel- 
eration with the same parameters. The higher Doppler-boosted parts of the flow may dominate 
the high energy emission of blazar jets. We suggest that highly variable GeV and TeV emission in 
blazars is produced by the faster moving leading edges of highly magnetized non-stationary ejec- 
tion blobs, while the radio data trace the slower-moving bulk flow. Model predictions compare 
favorably with the latest Fermi 7-ray and MOJAVE radio VLBI results. 



1. Introduction 
1.1. Bulk Lorentz Factor Crisis in AGNs 



One of the defining characteristics of many AGNs is their flux variability in all spectral bands ( Krolik 



1999). In particular, extremely short time scales of TeV variability are challenging to the models. The 



rapid flares reported for Mrk 501 and PKS 2155—304, on timescales of 3-5 minutes ( Albert|2007 Aharonian 



2007) imply an emitting size smaller than the gravitational radius tic ~hours of the supermassive black 
holes of these blazars. There are two contradictory issues related to short time scale variability. First, 
it implies a small emission size, which poses a problem for efficiency of energy conversion into radiation. 



Secondly, there is the compactness problem ( Guilbert et al.|1983 ). If variability is detected in 7-ray photons 



of energies exceeding the electron rest mass energy, then the emission region contains photons which can 
pair produce. If the number density of these photons is too high, then none of the photons will escape the 
region. The solution to both problems is bulk relativistic motion towards the observer, which reduces the 
intrinsic luminosity, decreases the implied energy of the photons, and increases the internal time scales. The 
required Doppler factor then exceeds 6 > 100. 

While highly relativistic motion may appear to be a cure-all, in AGNs the bulk Lorentz factor 7 can be 
directly constrained by VLBI observations of bright blobs moving with apparent speeds on the sky, /3app, that 
appear to be superluminal. This type of motion occurs when the emitting region is moving relativistically and 



close to the line of sight ( Rees 1966 ). The apparent transverse motion can exceed c due to propagation effects. 
If a blob is moving along with the bulk fiow of a jet and its velocity vector makes an angle, Oot, with the 
line of sight, then its apparent motion transverse to the line of sight will be: Papp — /3sin^ofc/(l ~ /3cos^ofc)- 
The maximum /3app can reach is /3"f when 9ob = I/7. Thus, if the blob motion corresponds to the underlying 
bulk motion of the jet, measuring /3app can constrain the possible bulk Lorentz factor, 7. 
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1.2. MOJAVE results 



The latest MOJAVE VLBI results do support the interpretation of moving jet features as physical 



entities, as opposed to patterns (Lister 2009a I. Observations of bidirectional motions, the near-absence of 



inward moving features, ejections of multiple blobs in the same jet with the same speed, and tight correlations 
of jet speeds with other properties, such as 7 ray emission, apparent 7-ray luminosity, brightness temperature, 
and even optical classification, all support the notion that the blob motion reflects the underlying flow. 

The MOJAVE survey of compact, highly beamed radio- loud AGN has analyzed the motion of emitting 
blobs in 127 jets and found that the observed superluminal speed distribution peaks at /3app ^ 10 and tapers 

off at jSapp 

around 



50 (Lister 2009b). This suggests that the bulk Lorentz factors of such objects are typically 
10, and extend up ^ 50, making the estimated values of 7 > 50 for PKS 2155—304 and Mrk 501 



rather difScult to reconcile with the radio data. Furthermore, direct VLBI observations of these sources on 
parsec scales have not even detected superluminal motion ( Finer &: Edward"s||2004 Giroletti et al.|[2004 |. 

VLBI observations of blazars such as PKS 2155—304 and Mrk 501 are not the only data which imply 
a low 7. Another way of investigating blazars is to search for their AGN counterparts whose jets are not 
directed along the line of sight, which are presumed to be radio galaxies (there are actually two distinct 
types of radio galaxies, FRI and FRII, that are thought to correspond to the two categories of blazars, BL 
Lacs and optically violent variable quasars Urry fc Padovani|1995 ) However, studies comparing the relative 
fluxes and numbers of radio galaxies and blazars point towards Lorentz factors of 7 < 10 (Henri & Sauge 
[2OO6). Indeed, preliminary results of MOJAVE observations of low luminosity jets from radio galaxies imply 
speeds of c or less. Thus, there is apparent contradiction between measured superluminal velocities and bulk 
Lorentz factors required by radiation modeling. This is known as the "Blazars' Bulk Lorentz Factor Crisis" 
dHenri fc Sauiil[2056 1 . 



2. Non-stationarity in AGN flows 



Following Blandford & Rees (19741, models of AGN jets typically assume steady-state injection con- 
ditions. This is based on the fact that the sonic time over the black hole horizon is typically shorter than 
any observed times scales of variability of black hole systems (except for such subtle effects as quasi-periodic 
oscillations in Galactic binaries, e.g., van der Klis 2005). By consequence, it is argued, a system reaches 



an quasi-equilibrium state. On the other hand, both Galactic black hole candidates, as well as AGN jets 
show a wide variety of non-stationary behavior. This non-stationarity is driven by various disk instabilities, 
occurring on the viscous time scale of (inner) accretion disk. 



The efficiency of the BH-powered jet depends {e.g., in a Blandford fc Znajek (19771 paradigm of jet 
launching) both on the parameters of the black hole (mass and spin) as well as on the magnetic field supplied 
by the disk. In addition, as we argue in this paper, jet acceleration may proceed more efficiently in case of 



non-stationary outflows (see also Lyutikov 2010, submitted, Granot et al.|[2010 ). 



One expects that accretion onto the black hole may change on the time scale of the order of viscous 
time scale of the inner part of the disk, which may be weeks to months to years for a typical AGN {e.g., 
if we associate blob ejection time scales as being due to the inner disk instability). As the relativistic jet 
propagates away from the central black hole it expands sideways, reaching a transverse scale of a fraction 
of parsec at the distance of the order of a parsec. Close to the black hole, the corona of the accretion disk 
is hot, with a nearly relativistic speed of sound: the sonic time scale of the magnetosphere across the black 
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hole is much shorter than the viscous time scale of the inner accretion disk. Further out, the variability time 
scale of the jet will be shorter than the dynamical time scale of the corona across the jet. Observationally, 



the connection of the disk and jet activities is well established in the case of Galactic microquasars (Mirabel 
fc Rodriguez||1999| ). In AGN, the data are controversial ( Chatterjee et al.]|2009 ). 



Let us assume that a jet with an opening angle of Qj ^ 0.1 is propagating through a corona, which is in 
hydrostatic equilibrium with the local sound speed Cg close to the virial velocity, = GMbh /'r- Relating 
the disk variability time scale to the Keplerian period near the black hole, — ^tbh/c ~ 10 days 9^2 j 
^ 3> 1, where the sound crossing time across the jet is longer than the variability time scale for 

/ f \2/3 

rbrcakout > J- j rsH = 2 X IQi^AfQ.gea/'eT^/^ cm. (1) 

where the notation Xa implies (X/IO") {e.g., ^2 = ■^/lOO, Oj,_i = 6^/0.1). At distances larger than Q, the 
jet variability on a time scale of a hundred Keplerian periods near the black hole, ^2 = 1 proceeds on time 
scales shorter than the dynamical time scale across the jet, so the external medium does not have time to 
react to changing jet conditions. It is also required that the jet travel time to rbrcakout be smaller than the 
the sound crossing time at rbrcakout- This requires ^ < l/0j- 

A newly created jet will then propagate along a nearly empty channel, cleared by the previous jet 
activity. Thus, we expect that initially, close to the black hole, the non-stationary injected jet propagates 
through a relaxed corona and will have to "bore" its way through. After reaching the distance ([T]), the 
leading edge of the jet will break out into low density medium created by the previous jet activity. What 
are the consequences of this non-stationary behavior for jet acceleration? 

As a model problem, we assume that a period of accretion lasting time has brought onto the black 
hole the magnetic field, initiating the Blandford-Znajek process. Since the external medium has fairly low 
density, the jet expansion is relativistic (see Eq. [9|. In a time td the jet inflates a bubble of linear size 
ctd ^ rbrcakout- At a distaucc rbrcakout, where the sonic time through the magnetosphere is of the order of 
the jet variability time scale, the newly inflated bubble reaches a near vacuum. What is the behavior of the 
jet starting at this point? 

As a model problem we considered (Lyutikov, submitted) a one-dimensional flow of cold magnetized 
into a vacuum and into an external medium of density pcx- This is reviewed in ^Js] In Qwe apply the results 
to non-stationary outflows in AGNs. 



3. Non-stationary relativistic expansion 

3.1. Riemann problem for relativistic expansion of magnetized gas 

Let us assume that before the breakout into the low density medium the jet plasma is moving with 
velocity (3^ towards the external medium. We have found (Lyutikov, submitted) an exact self-similar solution 
of the relativistic Riemann problem for the expansion of cold plasma with density po and magnetic field -Bq 
(magnetization parameter a — B"^! pQ\ the magnetic field is normalized by \/47r), moving initially with 
velocity towards the vacuum interface 
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6a 



.2/3cl/3 



(2) 



where the Doppler factors 5a = + /3a) /(I ^ /3a) are defined in terms of the plasma velocity /3, local 
Alfven velocity /3a, self-similar parameter i] = z/t, initial wind velocity /3m and the Alfven velocity in 
the undisturbed plasma /3a, o- These equations give the velocity /3, density p — U\pq/<7 (where Ua = 
(3a/ — P'a) and proper magnetic field, B — {p/po)Bo as a function of the self-similar variable rj — z/t 
(the expansion of plasma starts at t — 0, z = and proceeds into a positive direction z > 0). We stress that 
these solutions are exact, no assumptions about the value of the parameter a and velocity were made. 

Particularly simple relations are obtained for plasma initially at rest expanding into vacuum /3^ = 
0,5^ = 1 (Lyutikov, submitted). The flow accelerates from rest towards the vacuum interface (Fig. [ij). The 
bulk of the flow is moving with Lorentz factor 7' ~ a^^^. The flow becomes supersonic at rj = 0, at which 



point 7' = ((t/2)^/'^. The vacuum interface moves with Lorentz factor 7^,^^ 
the vacuum interface is moving with 5,, 



— I + 2a. In the observer frame 



6 A o^w, which in the limit cr, 7^ ^ 1 gives 



47m cr. 



(3) 



As the flow expands, the energy flux Too and the energy density Tqz stay nearly constant in the bulk of the 
flow at a value ~ Bq/A. The energy flux is maximal at the sonic point rj ~ 0, while the energy density slowly 
decreases towards the vacuum interface (Fig. [l]). 
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Fig. 1. — Energy density Tqq and energy flux Tq^ {Left panel) and Doppler factor = ^y {1 + f3) / {1 ~ (3) 
(Right panel) for self-similar expansion into a vacuum as functions of the self-similar coordinate 77 = z/t. 
Stationary initial conditions are assumed. The energy density is normalized to values in the undisturbed 
medium, while the energy flux is normalized to maximum values at rj — 0. Solutions extend from the front 

of the rarefaction wave at rjfiw — — •\/<7/(l + cr) up to the vacuum interface 77vac — 2 "^^{.2^)^'' • The Doppler 
factor is unity in the undisturbed plasma and reaches a value corresponding to 7^ac = 1 -I- 2cr at the vacuum 
interface. Plots are for a — 2. 

Most importantly, in a narrow region near r^vac = 2^/a{l + a)/{l + 2a) « 1 — l/(8cr^), the Doppler 
factor increases from ~ a^^^ in the bulk to Sp^max ~ 4(7 on the vacuum interface. Values of 6p > {l/2)Sp^max 
are reached within a range of A77 = 7/(8(T^) near the vacuum interface. The relative amount of energy 
with Doppler factors 6/} > (l/2)(5^,„iai is ^ S/cr^ (the total energy of the outflow at time t is Etot = 
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po{l + (t/2) ^J~j^t w Bgt/2). In addition, in a non-self-similar regime, most of plasma will reach the 
vacuum terminal velocity, since the forward characteristics never cross the vacuum interface regardless of 



the dimensionality of the flow (Greenspan & Butler 1962 Zeldovich & Raizer 2003 see also Lyutikov, 
submitted) . 



As the flow expands, the local magnetization 

P 



^2/3 

xU3 



decreases (Fig. [2]). At the sonic point aioc = (cr/2)^/^. 
10 F 




r2/3 
-'AM, 



(4) 




Fig. 2. — {Left panel) Local magnetization aioc, as given by Eq. Q. Initially a = 10. Right panel: Observed 
flux produced by the plasma expansion into vacuum, parametrized as a product of the rest frame energy 
density and the Doppler factor cubed. The scale is normalized to the value at 77 — 0. 



3.2. Expansion into medium 

If there is an outside medium with density pcx, we may identify two expansion regimes with different 
properties of the forward shock. For relativistically strong forward shocks, so that the post-shock pressure 
is much larger than density, the Lorentz factor of the contact discontinuity (CD) is 

(the last approximation assumes cr 3> 1). Here L Egj^c^r^ is the jet luminosity. For weak forward shocks, 
the velocity of the CD approaches the expansion velocity into vacuum 7„ac (Eq. [3| . The transition between 
the strong and weak shock occurs for 

g \ 1/3 

.204872 p, 

For a < (Jcrit , the forward shock is weak. In the limit cr 3> 1 the vacuum approximation is applicable for 

= ^ ^ 512.7>^c3.2 = 5 X 10 cm x (-) (-) , (7) 



^cr^t= [7^^^:^—] ■ (6) 
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where numerical estimates are given for typical parameters in AGNs. The required density is fairly low, yet 
it depends sensitively on the parameters of the flow. Even if the density is higher than given by Eq. ^, the 
resulting forward shock will propagate with a Lorentz factor that is only weakly dependent on the external 
density: 

JFS - 2^n,<j{Pex/Pe.fir^'^ = 200 X (^) (^) [Pex / Pe. fiT'^^ ■ (8) 



4. Implications for non-stationary acceleration in AGNs 

We propose that flares in TeV and GeV emitting blazars are produced in the leading expansion edges 
of non-stationary ejection events, moving with 7max ^ 100, while the observed velocities of the radio blobs 
correspond to the bulk motion with 7buik ^ 10 — 50. In this section we discuss observational implications of 
the model. 

In our picture, after the switch-on of the acceleration process, the AGN central engine produces a 
relativistic jet, and as the latter bores through the corona, it can reach a relativistic Lorentz factor 7^. The 
value of 7^, can be determined by the pressure balance (Eq. [5| with high external density: 

7. = f A) - 20L]%e,^^,n~'/\ (9) 

where we have assumed ^ = l/6j and typical parameters of AGN outflows. 

After breakout, at r > rbrcak, the total Lorentz factor in the observer frame will be 2a^/^ times larger 
for the bulk flow, ^ 2^yja^/^, and 4(7 times larger for the leading edge, ^ 47^,17: 

7,„i, = 40aV3iV4e,,_in-i/4 

7vac = 80 (7 4{fge,- (10) 

For values of a exceeding unity, the difference between Lorentz factors of the bulk flow and that of the 
leading edge is even greater. 

Soon after the break out, in the self-similar stage, the relative amount of energy in the fast leading tail 
is fairly small in the case of highly magnetized jets, ~ 3/cr^, yet it may dominate the beamed high energy 
emission due to high Doppler boosting. For example, if we parametrize the observed intensity produced by 
the jet as a product of the rest-frame energy density, Toz/ij^fi), and Doppler factor cubed (5|, it will be 
dominated by the fast moving parts of the flow (Fig. [2|. 

In addition, at later stages of expansion of a finite size blob, a larger fraction of the material may be 
accelerated to the maximal Lorentz factor. If the initial blob had a size L in its rest frame, most of the blob's 
material gets accelerated to Lorentz factor jI^^^^ ~ 1 -|- 2cr in time At ^ Sa'^L/c (Lyutikov, submitted). For 
a blob moving with 7^, the acceleration time in the observer frame is further extended to At ^ IGa^j^L/c. 



This time is, typically, longer than the time scale for blob interaction, Eq. (15). Thus, the bulk of the 
material does not get accelerated to the maximum Lorentz factor 7max- 

There are a number of correlations that we would expect in the model. First of all, we expect some 
correlation between the 7-ray and radio core fluxes, and this has already been seen by Fermi and MOJAVE 
( [Kovalev et aI1[2009l ) . 



Since the 7-rays are produced in a faster moving part of the flow, we expect that the jets of 7-ray selected 
AGNs are more aligned than those in radio-selected samples. This effect should be reflected in a flux-flux 
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plot of a well-defined flux-limited blazar sample, due to the different degrees of Doppler beaming in the radio 
versus 7-ray regimes. For example, a tight linear correlation of intrinsic (i.e., unbeamed) radio vs. 7-ray jet 
luminosity will be greatly smeared out in the flux-flux plot, although an upper envelope will still be present 
(see, e.g., simulations by Lister|[2007 ) . Indeed, the form and scatter of the observed radio-7-ray correlation 
in blazars ( Kovalev et al.|2009 ) does support the notion that the 7-ray emission is likely boosted by a higher 
Doppler factor than the radio emission. This is also reflected in the superluminal speeds and apparent jet 
opening angles, which suggest higher Doppler factors and smaller viewing angles for Fermi-detected AGN 
( [Lister et al.]|2009b[ [Pushkarev et al.||2009[ ) 



On the other hand, a radio-7-correlation depends crucially on the fact that the cores of AGN jets are 
optically thick to synchrotron emission up to the frequency-dependent radius ( Blandford fc Konigl|1979 Eq. 
28) 

-I . >2/3r2/3 -1/3 -1 -, n 

where (^j^ parametrizes the observed radio luminosity in terms of the total jet power Lj^ — CrL, with 
Cfl ^ 0.01; we also assumed here (and in the estimates below) that 8j ^ l/7ji C = The core radius 

( 11 1 is typically much larger than the breakout radius ([T]). Thus, typically, the jet breakout will occur while 
the jet is still optically thick in the radio. In this case the high energy emission will not be accompanied by 
the simultaneous increase of radio flux (at least at low frequencies). At sufficiently high radio frequencies, 
the radius where the jet becomes optically thin ( 11 ) may become comparable to the breakout radius ([!]). 



We expect that j-ray event should precede the radio blob ejection by 



At 



-f-R 



272 



,^2/3^2/3 7/3, 



days. 



(12) 



These are typical time scales between variability at radio, optical and 7-ray bands ( [Marscher et al. 20101 
Also, since this timescale (Eq. 12 ) is fairly long, we expect occasional 7-ray flares without radio blob ejection 



this appears to be the case in TeV flares associated with high-energy peaked AGN such as Mrk 421, and 
Mrk 501, in which no major changes in VLBI radio structure are seen (MO J AVE program; Ros et al. in 
prep). Furthermore, the ^-ray events should be better correlated with with radio blob ejection at high radio 
frequencies. Higher resolution monitoring of MOJAVE sources by the Boston University group (Jorstad 



et al.|[2005 ) at 43 GHz has revealed instances of jet features moving at higher speeds than those seen at 15 
GHz (Lister 2009a) that fade out very rapidly close to the core region. The extensive VLB A-|- Fermi data 



currently being gathered on these AGN should provide a useful test of this prediction. 

The above correlations assume that the 7-ray photons produced at the breakout radius ([T]) do not suffer 
from absorption. This depends on the compactness parameter corresponding to the bulk motion of the jet 
(since 7-ray photons can pair produce on a lower moving bulk plasma): 



I 



1 aT 



^e^^ ^breakout 



= lO^^L, 



7,440^ 



On the other hand, the 7-ray variability time scale depends on the Lorentz factor of the leading edge 

A , ^breakout 1 -i r-n ^ r\ 

At 2 9 = 150 sec Oj _i 

c 327jo-^ a'^ 



(13) 



(14) 



Thus, the model is able to accommodate both the requirement of small optical depth for j-ray photons and 
the short time-scale variably, down to few minutes. 
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The jets of lower power FRI sources and higher power FRII sources have somewhat different morphology 
on mas-scales. Relatively low power blazar jets in the MOJAVE program like Mrk 421 or Mrk 501 show 
a very smooth fall-off in radio intensity with distance from the core, and have maintained this structure 
over more than a decade of VLBI monitoring. In contrast, higher power blazars such as 3C 279 display 
jet morphologies dominated by individual bright knots that continually emerge from the core and move 
downstream at superluminal speeds (see Fig. [3|. 

One possible explanation for the dichotomy of jet properties is that the emitted blobs merge, creating 
smooth large scale profiles. If the blobs are ejected on a time scale td = (.tbh/c, they will merge at distance 
(in the observer frame the time of merger is determined by the Lorentz factor of the trailing edge of the 
preceding injection 7 '-^ 7j/ (4cr): 



2c*'^ll^)-3xl02"cm^e-!, (15) 



Thus, at linear scales > 100 pc, the jet is expected to be mostly smooth. Note, that since in blazars the 



angle between the jet direction and the line of sight are small, the projected distance corresponding to ( 15 ) 
are small. For example, with the jets in Mrk 421 or Mrk 501 at z=0.03 oriented at a few degrees to the line 
of sight, this corresponds to a few milliarcseconds projected on the sky. In addition, since the radius rmerge 
depends sensitively on the assumed bulk Lorentz factor, our model predicts that in powerful FR-II jets, the 
jet may remain knotty up to the kiloparsec scale. 

Can the observed morphology of the jets be used to determine intrinsic jet properties, like bulk Lorentz 
factor and magnetization? High-frequency peaked blazars (HBLs) are under-represented in the MOJAVE 



sample; in fact, there are no HBLs in the complete radio-selected sample (Lister & Homan 20051. On the 



other hand, HBLs are well represented in the Fermi sample fAbdo et al.|2009| In the framework of our model, 
this can be due to their lower bulk Lorentz factors, but high magnetization, see Eq. (10). It would then 
imply relatively smaller merging distances (15) and smoother jet structures, consistent with observations 
(Fig. 3). An obvious caveat in this argument is that HBLs have flatter 7-ray spectra and are more likely to 
be seen by the Fermi LAT. Further careful analysis of beaming and instrumental selection effects is needed 
to explore these possibilities more fully. 

FRI and FRII sources could also have different intrinsic variability times scales. In our present model, 
we relate the variability times scale to the mass of the central black hole, which does not vary considerably 
between FRI and FRII sources ([Ghisellini fc Celotti|[200T I . 



5. Discussion 

In this paper we discuss the effects associated with non-stationarity of the jet ejection. In particular, we 
argue that the leading edge of a non-stationary magnetized outflow can achieve a bulk Lorentz factor much 
larger than would be inferred for a steady state outflow given similar conditions. In the case of the expansion 
of a highly magnetized plasma, the ratio of Lorentz factors of the bulk flow and that of the leading edge 
can be as high as 2a^/^ {a is plasma magnetization). This ratio can reach tens for highly magnetized flows 
with a ~ 10. We suggest that the Doppler factor crisis in AGNs (a difference of the Doppler factors inferred 
from radiation modeling, especially of short time scale TeV flares and from the observations of radio blobs) 
may be resolved by non-stationary outflows: highly variable emission is produced by the fast-moving leading 
edge of an expansion, while the radio data trace the slower-moving bulk flow. 
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Fig. 3.— MOJAVE 15 GHz VLBA images of low-energy-peaked quasar 3C 279 (1253-055; top two panels), 
and high-energy-peaked BL Lac object Mrk 421 (11014-384; bottom two panels). The left-hand panels 
show the parsec-scale jet structure in a recent MOJAVE epoch. The right-hand panels show time-averaged 
images created by combining archival VLBA epochs from 1995 to 2009. The relatively smooth fall-off in jet 
intensity downstream from the bright core in Mrk 421 is a characteristic of high-energy-peaked jets in the 
MOJAVE sample. These jets display very few changes in their radio jet structure over time, as indicated by 
the similarity in the inner 5 milliarcseconds of the single epoch and stacked images of Mrk 421. By contrast, 
the jets of flat-spectrum radio quasars such as 3C 279 generally display a constantly varying structure that 
is dominated by numerous bright super luminal knots (Lister et al.||2009a Lister|[2009a). 



The suggested model is qualitatively different from the internal shock models, where non-stationarity 
is invoked to produce shocks and dissipate the energy of the relative bulk motion of the colliding media. 
Collision of strongly magnetized plasma blobs results in only weakly dissipative internal shocks. The fast 
leading expansion edges will generate powerful shocks in the surrounding medium that may produce the 
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high energy emission. We do not address the question of how magnetic and bulk energy is converted into 
radiation. 



A somewhat similar continuous acceleration mechanism was proposed by Tchekhovskoy et al. (2010) 



(see also Komissarov et al. 2009). It relies on sideways expansion of the jet after the break out. Sideways 
expansion of unconfined magnetically dominated plasma proceeds with Lorentz factor 1 + 2(T, so that the 
total Lorentz factor (of a plasma near the edge, affected by the rarefaction wave) is (1 + 2ct)7u,, two times 
smaller than for the case of expansion wave propagating along the direction of motion. This factor of 2 
may have an important effect on escape of high energy radiation, since the optical depth to pair production 
scales approximately as F"^, a difference in a factor of 2 in F will result in a difference of 64 in the optical 



depth. (Tchekhovskoy et al. (2010) cannot treat parallel acceleration akin to breaking into vacuum since the 



rotation of the central engine is turned on gradually for numerical reasons.) 

In summary, we believe that magnetically-driven non-stationary jet acceleration can provide a potential 
resolution of the longstanding bulk Lorentz factor crisis in blazars. The rich Fermi- VLBA dataset that is 
currently being gathered on a broad set of AGN should provide an excellent means of testing our proposed 
scenario. 
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